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Abstract. The deepest observations of the X-ray background approach the surface brightness of the truly diffuse 
component generated by Thomson scattering of cosmic X-ray photons. Available estimates of the electron density 
and the X-ray luminosity density of AGNs as a function of cosmological epoch are used to calculate the integral 
scattered X-ray background component. It is shown that the scattered component constitutes 1.0 — 1.7% of the 
total background, depending on the AGN cosmic evolution. Albeit this is a minute fragment of the total flux, it 
becomes a perceptible fraction of the still unresolved part of the background and should be taken into account in 
the future rigorous assessments of the X-ray background structure. This diffuse component at energies JS 1 keV 
sums up with the emission by WHIM to 3 — 4%. Consequently, one should expect that integrated counts of 
discrete sources account for just 96 — 97% for soft background and ~ 99 % at higher energies. 
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1. Introduction 



Most of the X-ray background ( XRB) is generated by 
discrete extragalac tic sources (e.g. ICampana et al. ifeoOlL 
iRosati et al. lEooil and references therein). Among those 
sources, various classes of AGNs constitute a dominating- 
part. Probably 5 to 10% of the soft XRB is produced 
by hot gas in clusters of galaxies. Hydrodynamical sim- 
ulations of the evolution of the primordial gas demon- 
strate that some contribution to the XRB is generated 
by the baryonic matter which fills intergalactic space in 
the form of moderately hot plasma llCen fc Ostrikerll 999. 
iDave et al. 1200(1 iBrvan fc Voit l200lllCroft et al. 120011) 
Due to its low density, thermal emission of this gas is 
extremely weak and the very existence of the "warm" 
bary ons still waits for d efinite observational confirmation 
(e.g. lSoltan et al. 2002). Thus, apart from this small con- 
stituent, it is generally accepted that the XRB is produced 
by discrete sources. This notion has been strengthened di- 
rectly by deep observations by RO 'SAT and more recently 
by Chandra. 

A relationship between the total XRB flux and in- 
tegrated flux produced by sources has be en investigated 
in det ail from the observational point by iMoretti et al. I 
(2003J. These authors compared the source counts with 
the measured XRB level in two energy bands: soft (0.5 — 
2keV) and hard (2 - lOkeV). They have shown that 
in the soft band the source counts integrated down to 
the Chandra limit of 2.4 x 10~ 17 ergs _1 cm~ 2 produce 
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94.3_ 6 ; 7 % of the XRB. A smooth extrapolation of counts 
down to ~ 3 x 10~ 17 ergs -1 cm~ 2 generates 96% of the 
XRB and is "consistent with its full value at 1 a level" . In 
the hard band discrete sources generate at most 93 % (af- 
ter source counts extrapolation) and are "only marginally 
consistent" with the XRB level. One should note also that 
th e total XRB flux itself is not well determined; according 
to IMoretti et al. I l)2003j) the uncertainties of the XRB in 
both bands amount roughly to 5%. 

Although the present data do not allow definite con- 
clusions, it is possible that the apparent deficiency of the 
source counts contribution to the hard band could be re- 
moved by steepening of counts below the present thresh- 
old of 2.6 x 10 _16 ergs" 1 cm~ 2 . Such steepening is in fact 
predicted b y recent evolutionary mo dels of the obscured 
AGNs fe.g. lFranceschini et al. Il2002|) . 

Since the deepest observations allow only for a few per- 
cent of the XRB in the form of genuine diffuse radiation, it 
is worthwhile to examine also those emission mechanisms 
which generate a weak signal as compared to the total 
XRB flux, but which potentially could contribute signifi- 
cantly to the diffuse component. The aim of the present 
paper is to estimate flux of the truly diffuse emission pro- 
duced by Thomson scattering of X-rays in intergalactic 
space. The amplitude of the integrated scattered flux de- 
pends on the cosmic history of the XRB and on the density 
of free electrons as a function of redshift. These functions 
are investigated in Sects. Inland respectively. In Sect.|H 
the amplitude of the scattered component is calculated 
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and it is discussed what constraints on the faint end of 
the source counts are imposed by the present results. 

Most of the AGN related XRB is generated at red- 
shifts smaller than ~ 3 with practically no contribution 
at redshifts greater than 6. Since cumulative Thomson 
scattering optical depths at redsh ift 3 and 6 do not ex- 
ceed respectively ~ 0.02 and 0.03 (|Cen II2003I) . the effects 
of X-ray absorption remain low and the calculations be- 
low do not include effects of multiple scattering of X-rays. 
Integration of the radiative transfer is reduced to calcula- 
tions of the XRB intensity as a function of redshift and 
then integration over redshift of the scattered component. 




2. AGN luminosity density and the X-ray 
background 

2.1. Local XRB 

The relationship between the average X-ray background 
intensity and redshift is investigated. The XRB flux in the 
local universe is know n with relatively high accuracy (e.g. 
iMoretti et al. 1 120031 and reference therein). It is estab- 
lished that a large fraction of the local XRB is produced 
by AGNs distributed over a wide range of redshifts. Due 
to the strong cosmic evolution of the AGN population, the 
XRB varied substantially in the past cosmological epochs 
and AGNs constituted a dominating source of the X-ray 
radiation for most of the lifetime of the universe. In the 
present estimates only the AGN contribution to the XRB 
is considered, although some other classes of objects such 
as young galaxies could contribute to the XRB at high 
redshifts. 

Observed locally at energy E the XRB flux I {E) is 
equal to the integrated flux produced by sources along the 
line of sight: 



h(E) 



dz 



dV C Z (E') ■ (1 + z) 



dz 



(1) 



where V is the co-moving volume, Dl is the luminosity 
distance, and C Z (E') is the luminosity density generated 
by AGNs at redshift z and energy E' = E ■ (1 + z). The 
maximum redshift, z max , is the redshift at which the first 
AGNs began to shine. Systematic variations of the lumi- 
nosity density with redshift are described by the cosmic 
evolution £(z): 



C z = C - £{z). 



(2) 



Various models for £(z) have been proposed in the lit- 
erature. The models differ between each other in details, 
but the general shape of the evolution is common. It is 
accepted that at low redshifts the luminosity density gen- 
erated by AGNs rises sharply, stabilizes at moderate red- 
shifts and probably decreases at high redshifts. One should 
note that in the present investigation the integrated lumi- 
nosity density as a function of redshift is the only "inter- 
esting" quantity. Thus, neither the exact shape of the X- 
ray luminosity function nor the evolution type (luminosity 



Fig. 1. Rate of the cosmological evolution of the luminos- 
ity density accord ing to the LDDE1 m odel with a decline 
given m Eq. (3) bvlMivaii et al. I l|200ri) (solid) and density 
evolution bv lBovle et al. 1 1 




(1 



(dashed). 



vs. densi ty evolution) a f fect th e calculations. It was found 
that the iMivaii et all l)2000|) estimates of the AGN lu- 
minosity function and cosmological evolution can be used 
to determine the luminosity density in a straightforward 
and effective way. Their Luminosity Dependent Density 
Evolution model (designated as LDDE1) is applied. It is 
based on the combined shallow and deep ROSAT surveys 
and adequately reproduces both the source counts and the 
AGN redshift distribution. 

The present calculations cover a wide range of red- 
shifts. The observed energy band [E\,E2\ corresponds to 
the band [E\(\ + z),E2(l + z)] at the source. Since the 
observed flux results from the integrated emission along 
the line of sight, the relevant luminosity function and evo- 
lution are d efined in the ( redshi fted) source frame (see 
discussion in lMivaii et al. I l)2000|) ). In the paper the soft 
X-ray band 0.5 — 2keV is used as reference and all the 
data extracted from the literature are corrected to this 
standard band. 

The AGN sample used in the IMivaii et al~l l|200f]h 
analysis contained a limited number of high redshift QSOs 
(17 with z > 2.2). In effect, constraints on the evolu- 
tion ra te at high redshifts are not restrictive. IMivaii et al. I 
(2000) note, however, that the X-ray data are consistent 
with the decline of the space density of optically selected 
QSO for redshifts greater than w 2.7. To account for this 
trend, the evolution function, £(z), is calculated from the 
LDDE1 model with an additional "damping factor" : 



£ z ~ exp(2.7 - z) for z > 2.7. 



(3) 



The cosmological model (f2 m ,f2A) = (0.3,0.7) has been 
used in all the calculations. In Fig. Q] the evolution rate 
defined as a ratio of the co-moving luminosity density at 
redshift z to its local value is shown for the LDDE1 mod el. 
The density evolution derived by Bovlc et al. is 
plotted for comparison. Differences between both evolu- 
tion laws of ~ 25% which arise at z 1.5 indicate the 
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level of uncertainties involved in the evolution estimates 
(see below). 

A relationship between volume, distance and redshift 
for the cosmo logical model used in the paper is taken from 
lHoggl<jl999|) : 



dV _L_ 

where 



H F Q {z){l + zf 



f {z) = ^n m {i + z) 3 + n A 



(4) 



(5) 



with fl m — 0.3 and £l\ = 0.7. Eqs. (4) and (5) inserted 
into Eq. (1) give for unit solid angle u> = 1 sr: 



I (E) = -L-^£ (E) 



dz 



£{z) 



F (z)(l + z) 2 



(6) 



The local background flux I (E) is only weakly de- 
pendent on the maximum redshift in the integral and 
in all the calculations z max = 6 was put. The lo- 
cal luminosity density and e v olutio n have been com- 
puted from the iMivaii et al. I l|200fj|) formulae and in- 
serted into Eq. (6). The results are dependent on ex- 
trapolation of the X-ray luminosity function (XLF) at 
the low luminosity end. In the present calculations the 
XLF is integr ated over the l umino sity range actually cov- 
ered by the IMivaii et al. I l)2000|) sample, i.e. between 
Lmin = 10 41 ' 5 /i^ 2 ergs" 1 and L max = 10 47 ^ 2 ergs _1 , 
where /150 = iJ o /50kms _1 Mpc~ 1 . In the energy band 
of 0.5 - 2keV the AGNs according to the LDDE1 with 
the decline described by Eq. (3) model produce a back- 
ground flux of 1.97 x 10 _8 ergs~ 1 cm" 2 sr" 1 . Using the 
iMoretti et al~l l)2003|) estimate of the total soft XRB at 
2.47 x 10 -8 ergs -1 cm -2 sr -1 , it gives the AGN contribu- 
tion at ~ 80%. A pure LDDE1 model (without decline of 
luminosity density at high redshifts) generates a flux of 
2.16 x 10- 8 ergs" 1 cm" 2 sr" 1 (= 88% of the total XRB). 



2.2. XRB in the past 

A hypothetical observer at redshift z would detect the 
integrated XRB appropriate for his epoch: 



I'AE' 



, dV> C' Zl [E'(l + z')]-(l + z') 

dz 17 ^ ' (7) 



4nD' 2 



where primes (') indicate that all quantities refer to and 
are measured by the observer at redshift z. Photons emit- 
ted at redshift z' were sent at an epoch, which for us cor- 
responds to the redshift z\\ 



1 + *i = (l+ *)•(! + *') 



(8) 



The cosmological parameters in Eq. (7) satisfy relation- 
ships analogous to Eqs. (4) and (5): 



dV 1 c 
u> 



dz' D 2 / H z F z {z'){l + z>) 



A 2 



(9) 




Fig. 2. XRB flux in the band 0.5(1 + z)keV - 2.0(1 + z)keV 
per co- moving cm 2 measured at redshift z. The solid curve 
represents data for the LDDE1 model with decreasing lu- 
minosity density according to Eq. (3), the dotted curve - 
the LDDE1 model without "damping". In both cases the 
luminosity density is set to zero for z > 6. 



and 



F z (z') = yJn> m (i + z>)s + n' A . 

Standard relationships (e.g. iHogg lll999]) give: 
H z = Ho-F {z) , 
F ( Zl ) = F (z) ■ F z (z') , 
£' zi = (1 + zfC zl 
1 



h(E') = 



(1 



■I'AE'), 



(10) 

(11) 

(12) 
(13) 

(14) 



where I' Z (E') is the XRB flux per proper unit surface mea- 
sured by observer at redshift z and I Z (E') represents flux 
per co-moving unit surface. Using Eqs. (8) - (14) we finally 
get: 



J Z (E[,E' 2/ 



— (l + z)Co(E 1 ,E 2 ) 
-Ho 



dz' 



F ( Zl )(l + z') 2 ' 



(15) 



where J Z (E[, E' 2 ) is the XRB flux per co- moving unit sur- 
face measured at redshift z within the energy band E'^—E^ 
multiplied by a solid angle of 4-7T. 

3. Electron density in intergalactic space 

Extensive hy drodynamical simulat ions of the matter dis- 
tribution fe.g lValageas et al. 120 02 and references therein) 
show that a dominating fraction of all the baryonic mat- 
ter in the universe resides outside galaxies and cluster 
of galaxies. In the local universe baryons accumulated in 
these objects constitute ~ 30 % of the total baryonic mass, 
while the remaining ~ 70 % still stays in the diffuse and 
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warm-hot phases filling 99 % of volume space. The fraction 
of the diffuse compone nts rises with redshift and reaches 
rough ly 99% at z = 3 l|Cen fc Ostrikedll999l iDave et al. I 
2000). Inte r galactic m atter is virtually fully ionized up to 

z «fi fR.g.El bnnah . 

Estimates of the total baryonic mass based on WMAP 
observations give Slh = 0.044 for H — 71 kms" 1 Mpc -1 
ijSpergel et al. 1 120031) . Using results on the fraction of 
baryons in the intergalactic space by ICen fc Ostrikerl 
l)l999t) . the WMAP data give the local average electron 
density at n e (0) = 2.2 x 10 -7 cm -3 . Due to systematic 
increase of the diffuse component with redshift, the elec- 
tron co- moving density rises to n e (z) = 2.2 x 10 -6 cm -3 
at z = 1 and to n e (z) = 2.0 x 10 -5 cm -3 at z = 3. At 
still higher redshifts, practically all baryons reside in the 
diffuse component and n e (z) ~ (1 + z) 3 . 



4. The scattered component 

The luminosity density generated by the Thomson scat- 
tering of the XRB is given by: 



Cl = n e {z) o~t J z 



(16) 



where ctt = 6.65 x 10 -25 cm 2 is the Thomson cross-section. 
The XRB flux due to Thomson scattering is calculated 
using an equation analogous to Eq. (6): 



^{E U E 2 ) = ^±- 
Air Ho 



dz 



F (z)(i + Z y 



(17) 



Due to the strong dependence of the electron density 
on redshift, the scattered component is generated primar- 
ily at high redshifts. The standard LDDE1 model with 
constant luminosity density for 1.6 < z < 6 represented 
by the dotted curve in Fig. [2] generates scattered XRB 
flux of 4.14 x 10 -10 ergs -1 cm -2 sr -1 , which accounts for 
1.7% of the total XRB. If the AGN luminosity density 
declines exponentially for redshifts z > 2.7 according to 
Eq. (3) (solid curve in Fig. i the scattered XRB flux is 
reduced to 2.42 x 10 -10 ergs -1 cm" 2 sr -1 or 0.98% of the 
total XRB. 

Low optical depth for Thomson scattering for z < 6 
implies relatively weak signal of the scattered XRB com- 
ponent. However, in the deepest exposures of the present- 
day instruments just a few percent of the soft XRB re- 
main unresolved into discrete sources. Thus, the scat- 
tered component contributes measurably to the unre- 
solved part. Moreover, a separate diffuse contribution by 
the WHIM is expected in the soft XRB. The amplitude 
of this component has not been determ ined yet with rea- 
sonable accuracy l|Soltan et al. Il2002f) . The WHIM con- 
tribution to_th£_JCRB_is expected to drop sharply with 
redshift l|Cen fc Ostrikerl lT999) and it is likely that ther- 
mal emission in the 0.5 — 2.0 keV band does n ot exceed 
2 - 3 % of the total XRB. IMoretti et al. I l)2003ft point out 
that a smooth extrapolation of the counts down to a flux 
level of ~ 3 x 10 -18 ergs -1 cm -2 sr -1 , i.e. a factor of 10 
below the present limits could account for 96 % of the to- 
tal XRB. Although this estimate is subject to relatively 



high uncertainties (~ 7% at the la level), the existence 
of the diffuse component further strengthens the conclu- 
sion that with a moderate extrapolation of source counts, 
all constituents of the XRB are accounted for and there is 
no room for new classes of X-ray sources. In particular, in 
the soft band the source counts cannot exhibit steepening 
which would increase the integrated XRB more than a few 
percent. 

In the harder band of 2 — lOkeV source counts ex- 
trapolated down to fluxes an order of magnitude be- 
low t he present threshold generate ~ 93 % of the total 
XRB (|Moretti et al. Il2003|) . The relative contribution of 
the Thomson scattered component in this band is simi- 
lar to the contribution in soft X-rays, while the expected 
WHIM emission above 2keV is negligible. Thus, taking 
into account t h e diffu se component does not change the 
IMoretti et al. ] ll2003h conclusion that the source counts 
are likely to steepen below ~ 10 -15 ergs -1 cm -2 due to a 
population of highly obscured, hard sources. 

Taking into account the uncertainties of the absolute 
flux of the total XRB and uncertainties of the discrete 
source contribution, which both are at the level of several 
percent, calculations in the present paper do not address 
the question whether the source counts match exactly 
the observed XRB. One should stress that the objective 
of the present investigation was to estimate the amount 
of the diffuses XRB component generated by Thomson 
scattering. The amplitude of this component is calcu- 
lated using the AGN X-ray luminosity function and evo- 
lution. The uncertainty of the diffuse component flux re- 
sults mainly from uncertainties related to the evolution 
of the AGN population. The present calculations show 
that the Thomson scattered flux amounts to or slightly 
exceeds 1 % of the total XRB, Thus, as long as uncertain- 
ties of the relevant measurements are subject to larger 
errors, this component does not affect decisively the XRB 
budget. Nevertheless, in the soft band the Thomson com- 
ponent narrows noticeably the distance between the re- 
solved and total XRB and reduces the range of the al- 
lowed logS* — logA^ relatioships. One might expect that 
the relative importance of the Thomson component will 
grow with a further increase of the resolved fraction of the 
background. 

Direct detection of the scattered component poses a 
serious observational problem. This is because the spec- 
tral shape of the scattered flux strictly corresponds to the 
average AGN spectrum and the diffuse component mimics 
a population of unresolved AGN-like sources. The ampli- 
tude of the diffuse flux could be determined only indirectly 
by subtraction of the discrete source contribution from the 
total XRB. In practical terms, to establish precisely the 
flux produced by discrete sources one needs to determine 
the log N — log S relationship well below the present limits. 
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